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ABSTRACT 

There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by 
fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous 
plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal 
ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse 
bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the under- 
lying velocity distributions of the interacting species. In this paper, we identify the conditions under 
which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in 
shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence 
favors electron temperatures almost comparable to, and certainly not very deficient relative to proton 
temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best 
a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in 
most models invoking shock acceleration in discrete sources such as supernova remnants. However, on 
scales £ 100 pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray elec- 
trons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly 
detectable. 

Subject headings: acceleration of particles — cosmic rays — supernova remnants — radiation 
mechanisms: non-thermal — gamma-rays: theory 

1. INTRODUCTION 

The process of inverse bremsstrahlung has received attention from time to time over the last three decades, in various astrophysical 
settings. Inverse bremsstrahlung is defined here to be the emission of a single photon when a high speed ion collides with an electron 
that is effectively at rest, and has often been referred to as suprathermal proton bremsstrahlung (e.g. Boldt k Serlemitsos 1969; 
Brown 1970; Jones 1971; Haug 1972). Being a kinematic inverse (an analogy can be drawn with inverse Compton scattering) arising 
purely via Lorentz transformations between reference frames (i.e. it is still the electron that radiates), it is not to be confused with 
the true quantum electrodynamical inverse of conventional bremsstrahlung, namely the absorption of photons in 3-body collisions 
with electrons and ions, a process that is highly improbable in tenuous (i.e. optically thin) astrophysical plasmas. 

The first application of suprathermal proton bremsstrahlung in astrophysical models dates to the work of Hayakawa k Matsuoka 
(1964), which considered such radiation in collisions between cosmic rays and ambient electrons in the intergalactic medium as a 
source of the cosmic X-ray background (XRB) at energies ^ lOkeV. This concept was developed in the papers by Hayakawa (1969), 
Boldt k Serlemitsos (1969), and Brown (1970), with Boldt k Serlemitsos proposing that suprathermal proton bremsstrahlung could 
account for the flat spectral index of the XRB below around 15 keV, while noting that an unusually high column density of ambient 
electrons would be required to match the observed flux. The interest in the relevance of inverse bremsstrahlung to the cosmic XRB 
has since diminished, and more recent opinion favours the accumulation of unresolved discrete sources as the dominant contributions 
(see Fabian k Barcons 1992, for a comprehensive review). Although thermal bremsstrahlung at around 40 keV can fit the 2-100 
keV spectrum fairly well (Marshall, et al. 1980; see also Boldt 1987, for a review), recent critical evidence limits the contribution of 
truly diffuse emission from the intergalactic medium to less than 0.01% of the observed flux. This bound is derived from upper limits 
to the mean cosmological electron density obtained from the measurements of COBE to Compton scattering-induced distortion of 
the cosmic microwave background (Wright, et al. 1994). 

Renewed interest in inverse bremsstrahlung has arisen in the last year. Tatischeff, Ramaty k Kozlovsky (1998) and Dogiel et al. 
(1998) have discussed expectations for suprathermal proton bremsstrahlung in the X-ray band, assuming a significant enhancement 
of the density of low energy cosmic ray ions in the Orion region. Inferences of such an abundance of ions were drawn from the 
reported (e.g. Bloemen et al. 1994), but recently retracted (Bloemen et al. 1999), detection with COMPTEL/CGRO of emission 
lines from the Orion molecular cloud complex; these were attributed to Carbon and Oxygen nuclear line emission in collisions with 
ambient ions. In a different context, Valinia and Marshall (1998) conjectured that inverse bremsstrahlung from cosmic rays could 
be responsible for much of the diffuse X-ray emission seen in the galactic ridge. Pohl (1998) has responded to this proposition by 
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arguing that it is difficult to produce the diffuse component in the ridge via inverse bremsstrahlung without violating observational 
limits to nuclear excitation line and pion decay continuum emission. Furthermore, Tatischeff, Ramaty & Valinia (1999) argue 
that by normalizing the low energy cosmic ray ion populations using the galactic Be production rate, the contribution of inverse 
bremsstrahlung to the thin disk component of unresolved/diffuse 10-60 keV emission detected by RXTE is at most a few percent . 

The focus of the presentation here is on collisionless shocked astrophysical plasmas, which distinguish themselves from the 
diffuse emission scenarios addressed by Valinia and Marshall (1998). Pohl (1998) and Tatischeff, Ramaty & Valinia (1999) by 
the insignificance of Coulomb collisions compared to interactions between the magnetic field and charged particles. The principal 
examples in mind here are young supernova remnants, i.e. those perhaps mature enough to be in the Sedov epoch, but well before 
the radiative phase. The bremsstrahlung emission properties of shocked environments are dictated by the dissipation between ions 
and electrons in the shock layer. We show in Section 2 that only when electron and ion thermal speeds are comparable (i.e., when the 
temperature ratio is on the order of, or less than the mass ratio, T e /T p ^ m e /m p ) can inverse bremsstrahlung contribute significantly 
relative to normal electron bremsstrahlung in the optical, X-ray and gamma-ray wavebands. We argue in Subsection 2.4 that such 
conditions are not easily realized in the majority of shocked astrophysical environments, thereby rendering suprathermal proton 
bremsstrahlung unimportant for most discrete sources. By extension, the importance of inverse bremsstrahlung as a continuum 
emission mechanism in extended regions is contingent upon there being a relative paucity of primary cosmic ray electrons in the 
interstellar medium on appropriate length scales (^ 100 pc). This point, and the role of knock-on electron bremsstrahlung, are 
discussed in Section 3. 


2. INVERSE BREMSSTRAHLUNG VERSUS ELECTRON BREMSSTRAHLUNG 

The key question of interest to astrophysicists that we address here is when the inverse bremsstrahlung process is significant 
relative to the normal bremsstrahlung mechanism in shocked astrophysical plasmas. Clearly, the relative numbers of projectile 
and target particles will be central to answering this question. Here we aim to determine how such numbers and the associated 
emissivities for the two processes depend on standard plasma parameters. The emissivities, or photon production rates for either 
process, can be written in the form 

— = n,jf dpn CR {p)c0^j , (1) 

for shock-accelerated cosmic ray (projectile) particles with momentum distribution n CR (p) colliding with targets of density n* that 
are effectively at rest in the observer’s frame. Here do/du is the cross-section, differential in the photon energy, where we adopt 
the convention of using dimensionless units throughout this paper: photon energies E ^ are expressed in units of the electron rest 
mass m e c 2 , i.e., via u = E 1 /(m e c 2 ) . Hence do/du has c.g.s. units of cm 2 . Also, c/? is the speed of the projectile of momentum 

p . 

The distribution functions of accelerated particles in shock environs need to be described numerically, in general, as will be 
evident from the Monte Carlo-generated distributions illustrated later in the paper. However, we can provide a good estimate of 
the relative importance of inverse bremsstrahlung and classical bremsstrahlung using simple but representative analytic forms for 
the particle distribution functions. For these purposes, let us approximate the cosmic ray (shock- accelerated) proton and electron 
distributions in a discrete source by power-laws in momentum, broken at thermal energies. Defining cfl TiP and c/? Tie to be the 
proton and electron thermal speeds, respectively, their momentum distributions can be cast in the approximate form 


n »{p) = 


Pr,i 


3(r. - 13- 

r, - 1 + 3e, 



P < Pt.j 
P > Pt,j 


( 2 ) 


where p T>3 = m s cf3 T}S is the thermal (non-relativistic) momentum for species s (= e,p). Here, n 3 represents the total density and 
T s is the power-law index of a given projectile species. This discontinuous “schematic” form is chosen to mimic the slopes of the low 
momentum portion of a Maxwell-Boltzmann distribution and the non-thermal suprathermal portion of the cosmic ray distribution, 
and hence does not accurately describe the spectral structure near the thermal peak. Consideration of such spectral details, which 
depend on the nature of dissipation in the shock layer, would introduce only modest corrections to the rates derived below, and are 
immaterial to the qualitative conclusions of this paper. The factor e s is introduced to account for the fact that the non-thermal 
tail does not extrapolate directly from the thermal peak in shock acceleration-generated populations, and its values £ I define 
the efficiency of the acceleration mechanism, which can be a strong function of the field obliquity (e.g. Baring, Ellison t Jones 
1993) and the strength of particle scattering (e.g. Ellison, Baring & Jones 1996) in the shock environs. For ions, acceleration is 
generally quite efficient, implying e p ~ 0.1 . For electrons, the situation is less clear, pertaining to the well-known electron injection 
problem in non-relativistic shocks. However, as will be mentioned in Subsection 2.4 below, the observed electron-to-proton cosmic 
ray abundance ratio provides a canonical lower bound to the injection efficiency of electrons, and generally implies that e e cannot 
be dramatically less than e p . Finally, note that the target density of cool particles for the two processes in question will be written 
as Tit > and will represent a subset of either of the accelerated populations that are formed from the shock heating of the interstellar 
medium. 


2.1. The Bremsstrahlung Ratio in a Nutshell 

The ratio 7£ 0 -X of inverse bremsstrahlung to bremsstrahlung emission in the optical to X-ray and soft gamma-ray bands can 
be quickly written down without the encumbrance of the mathematical complexity of the differential cross-sections involved. Here 
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we enunciate this, the principal result of the paper in a simple and enlightening manner, so that the detailed expositions of the 
next two subsections can be bypassed by the reader, if desired. At energies below around a few hundred keV, the waveband of 
interest to the discussions of Boldt & Serlemitsos (1969), Hayakawa (1969), Brown (1970), Tatischeff, Ramaty & Kozlovsky (1998), 
Dogiel et al. (1998) and Valinia and Marshall (1998), the non-relativistic differential cross-section in Eq. (4) below is operable, and 
applies to both bremsstrahlung and inverse bremsstrahlung. This cross-section is a function only of the speed of the ballistic particle 
involved, for a given photon energy. Hence, any disparity in the emissivities for the two processes can only be due to differences 
in the numbers of ballistic protons and electrons at a given speed, a property that does not extend to relativistic domains. Since 
bremsstrahlung X-rays in supernova remnant shocks are generally produced by slightly superthermal electrons, this number is just 
n e Ce/P Te , reflecting the efficiency of injection into the shock acceleration process. For protons, whose thermal speeds are much 
lower, only a small fraction can participate in inverse bremsstrahlung collisions with shocked electrons, namely those with ^speeds 
exceeding the electron thermal speed /3 T>e . Using the suprathermal portion of Eq. (2), this constitutes a “density” n p c p /3^ p /? T ,e P , 
so that it immediately follows that the ratio of the two processes is of the order of 


R-o-x 



(3) 


for n e ~ n p , as required by charge neutrality. Therein lies the principal result, borne out in the derivations of the next two 
subsections, which are expounded because of their usefulness in astrophysical problems. Clearly, the ratio Ho-x depends principally 
on the degree of temperature equilibration, or otherwise, in the shocked plasma. 

2.2. Radiation from Non-relativistic Particles 

A more-detailed estimate of the relative importance of inverse bremsstrahlung and classical bremsstrahlung at optical to X-ray 
energies can be obtained by considering only non-relativistic accelerated ions, i.e. those with (5 1 . In this regime, the differential 

cross-section da/du for either bremsstrahlung or its inverse can be obtained from textbooks such as Jauch and Rohrlich (1980), 
and is the non-relativistic specialization of the Bethe-Heitler formula: 


da 

cfw 


= ¥-z 



2ira f Z ^ < 1 


(4) 


where Of = e 2 /(hc) is the fine structure constant, +Ze is the nuclear charge, and ro — e 2 /(m e c 2 ) is the classical election radius. 
Here is the relative speed (in units of c) between the projectile and target, namely the electron speed for classical bremsstrahlung, 
and the proton (or ion) speed for the inverse process. The applicability of this formula, which is integrated over final electron and 
photon angles, for both bremsstrahlung and its inverse follows from the Lorentz invariance of the total cross-section a and the 
effective invariance of photon angles and energies in non-relativistic transformations between the proton and electron rest frames. 
With this interpretation, c\J fi 2 - 2w represents the final electron speed in the proton rest frame for either process. Note that Eq. (4) 
is derived in the Bom or plane-wave approximation, and is strictly not applicable to regimes where /? £ 27rafZ or u> (3 <1 
where Coulomb perturbations to the electron wave-functions become important (see the end of this section). Appropriate results 
for low frequencies and classical regimes can be found in Gould (1970), and references therein, and amount to modifications of the 

argument of the logarithm. ...... 

Using Equation (1), it is straightforward to write down production rates of photons for the two processes m the limit ot non- 
relativistic projectile speeds. For regular bremsstrahlung, the vast majority of shock- accelerated protons are effectively cold due to 
their generally low thermal speed (except for unusually hot proton components in tenuous plasmas), so that the proton density n p 
represents the target density n t . Hence the photon production rate is 


dn-y(u) 

dt 


ibrems 


16 n p Tl e ^ 

T Pr,e 


c 

w 



r c ) 4“ Ce fi T.e, 



(5) 


where T e is the electron power-law spectral index, and 


= £{(/?-» 




( 6 ) 


(3 C = max | /? r , . (n 

Here we use the argument T to represent either T e or T p , as the case may be, in the work below. The first integral in Equation (6) 
specializes to [log e (2 #/w) - l]/2 in the limit of » u . The second integral in Equation (6) is generally expressible in terms 
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of hypergeometric functions, though such a manipulation is not enlightening. Specialization to the 0 2 w limit yields tractable 
integrals and the result [log c (2 /??/<*>) T 2/T]/T . However, remembering that modest corrections to the cross-section in Equation (4) 
are required (Gould 1970) in this low frequency limit, we restrict use of this specialization to the range 0\ £ . The other interesting 

limit is when 1 , for which only the second integral contributes; the result in this case is ^(w, 0 T , T) « 2 r “ 2 B(T , T) , 

where £(x,z) is the beta function. The substitutions 0 = \J%j cosh0, and identity 3.512.1 in Gradshteyn and Ryzhik (1980) 
facilitate these developments. The limiting cases can then be summarized as 


c/n 7 (w) 

dt 


brems 


a r e - 

0T,e E e — 1 + 3c e 



« £ Pl.e « 1 . 


( 8 ) 


where T(x) is the Gamma function, and 

A = Z 2 n p n e afr\c (9) 

defines the fundamental scale for the bremsstrahlung emissivity (i.e. at u> ^ 1 and 0 T — 1). The w” 1 behaviour at low energies 
reflects the convolution of the bremsstrahlung infra-red divergence with a distribution possessing a finite number of particles. The 
^-(i+r^/2) dependence above electron thermal energies arises because the bremsstrahlung flux spectrum traces the electron energy 
distribution. Note that Eq. (8) is strictly valid only when T e > 1 , which will always be the case for convergent electron distributions. 
A simplified derivation of the form of Eq. (8) is presented in the Appendix. 

The inverse bremsstrahlung spectrum can be evaluated in a similar manner, but noting that zero contribution arises whenever 
the proton speed drops below the electron thermal speed; such a domain of phase space corresponds to normal bremsstrahlung. 
Hence, thermal protons are never sampled unless the proton thermal speed is comparable to that of the electrons. It follows that 
the emissivity resembles the form in Equation (5) with just an / 2 -type term; extracting the appropriate power-law factor, the 
inverse bremsstrahlung emissivity 7 can be written 


dn 7 ( up 


dt 


inv 

brems 


16 n p n e ^2 &f r oC 3(r p 1) / 0t,p\ f r 0 p \ 


( 10 ) 


where T p is the proton power-law index. The limiting forms of the / 2 function are just as for the bremsstrahlung case, so that we 
quickly arrive at the limits 


dn^{Lj) 

dt 


16 


inv 

brems 


0T,e Tp IT 3e p /7x,e ^ 




rw 


2r ' /2 Ji fat , Pl.t < W « 1. 


( 11 ) 


c p i+r 
ur T 1 


The two photon ranges of interest here, namely lj & 0 2 e 1 and 0 2 e < 1, are the same as those for bremsstrahlung. This 

is due in part to the fact that only protons with speeds above electron thermal speeds will participate in inverse bremsstrahlung 
interactions. Note that relaxing the 0 P /? TjC requirement to bounds on j3 p greater than 0 TtC would introduce the appropriate 
numerical factor in Equation (11). Clearly, the factor [e p /e fi ](0 ] t#/ P r } e) Tp ~ l roughly defines the ratio of inverse bremsstrahlung 
to bremsstrahlung emissivities; expectations for its value in shocked plasmas will form the center of the discussion in Section 2.4 
below. 

It is salient to remark at this point that the differential cross-section in Equation (4) that we are using (and also the ultra- 
relativistic limit in Eq. (14) below 7 of the Bethe-Heitler result) was obtained from QED calculations in the Born approximation. At 
non-relativistic speeds, the Coulomb potential of the target electron or ion perturbs the projectile electron wave-function sufficiently 
that the plane-wave approximation breaks down. This occurs when p € x drops below % , where p e = 0m e c is the momentum of 
the energetic electron. Here x is the spatial scale appropriate to the interaction, which is the classical electron radius ro . Hence, 
so-called Coulomb corrections become necessary when 0 & ocf = e 2 /(he ) . In such cases, the matrix element for scattering must 
be determined using Coulomb wave functions. This has been done exactly for non-relativistic electrons colliding with nuclei by 
Sommerfeld (1931), leading to the application of a simple corrective multiplicative factor (Elwert 1939) 


Cse{0, w) = 4r 


1 — exp 

—2notiZ/0 

1 — exp 

-2*a,Z/p' 


, ff = y/0*-2u> , 


known as the Sommerfeld- Elwert factor, to the Bethe-Heitler cross-section: 
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Here is the electron’s speed (in units of c) in the nuclear rest frame after collision, and Z is the charge number of the nucleus. 
Note that a similar factor can be employed for electron-electron collisions (Maxon and Corman 1967, see also Haug 1975). 

Such factors become important for projectile electron speeds below around c/10 (for proton targets). When /? T|C < 2tv a r Z , 
the Coulomb correction factor is simply 0/y/0 2 - 2u> , and it propagates through all the integrands of the above developments. It 
is then quickly ascertained that in the w £ /?? >c < 1 limiting cases, the correction factor is close to unity and therefore can be 
neglected (modest low frequency corrections for w /? 2 are discussed in Gould 1970). For higher photon energies, the Coulomb 
corrections should be included, with the appropriate modification to the algebraic manipulations. While the resulting integral is 
not tractable except in terms of higher order hypergeometric functions, we determined numerically that the corrective factor that 
should be applied to the e « w « 1 cases in Equations ( 8 ) and (11) increases monotonically from 1.17 to 1.94 as T e or T p 
vary between 1 to 5, the representative range of power-law distribution indices. Since such factors are of the order of unity, we can 
assert that Coulomb correction factors are only marginally important for bremsstrahlung computations, and will prove immaterial 
to the considerations and conclusions of this paper; we neglect them in the algebraic expressions hereafter. 


2.3. Emission from Ultra-relativistic Species 

Hard X-rays and gamma-rays from bremsstrahlung and inverse bremsstrahlung are generated by ultra-relativistic electrons and 
ions. In this regime, a single cross-section cannot be used for the two mechanisms due to the introduction of aberrations in photon 
angles and modifications of energies in transforming between the two rest frames of the interacting particles. For normal e - p 
bremsstrahlung, the differential cross-section is the ultra-relativistic specialization of the Bethe-Heitler formula (e.g. see Jauch and 
Rohrlich 1980), obtained in the Born approximation: 


do 


du> 




27(7 -up 
u? 



(14) 


and is not subject to significant Coulomb corrections. Here, 7 is the electron Lorentz factor. This formula can be used to compute 
contributions to the bremsstrahlung spectrum from mildly relativistic energies upwards. The review paper of Blumenthal and 
Gould (1970) provides a useful and enlightening discussion of various features and issues of bremsstrahlung from ultra-relativistic 
particles. From Equation (2), assuming a non-relativistic thermal speed c/? Ti , for a species s, the Lorentz factor distribution is 


easily deduced to be 


My) 


Ksfs 3(r 5 - 1) _r. 

(A-.,) 1 -*'* IW + 3e, 7 


(15) 


This can be readily folded with the differential cross-section in Equation (14) via an adaption of Equation (1) to yield an emissivity 
appropriate to the gamma-ray band: 
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where A is given by equation (9), and 
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where rp(x) = c^log*, T(x)]/dx is the logarithmic derivative of the Gamma function, and C = ip{l) = 0.577215 ... is Euler’s constant. 
The evaluation of the integral is facilitated by identity 4.253.6 of Gradshteyn and Ryzhik (1980). The u ** 1 values of Equations ( 8 ) 
and (16) are dearly of comparable magnitude. Note also that the bremsstrahlung emissivity traces the electron energy distribution 
in the relativistic limit. A simple derivation of the form of Eq. (16) is provided in the Appendix. 

Determining the emissivity for inverse bremsstrahlung in the limit of ultra-relativistic ions is, in principal, more involved, largely 
because there are no published expressions for the differential cross-section, integrated over the various interaction angles, that 
are general enough to take the place of the Bethe-Heitler formula. Exact numerical computations in the Born approximation 
are presented by Haug (1972). Simple analytic formulae for the cross-section were obtained by Jones (1971), who used the well- 
known Weizsacker-Williams method (Weizsacker 1934; Williams 1935) where the process is treated as stationary electrons Compton 
scattering the virtual photons carried by the proton’s electromagnetic field. Such a technique is quite applicable to hard X-ray and 
soft gamma-ray energies, but becomes erroneous when u; > 1 since then the approximation of the Coulomb potential of the proton 
by a Fourier decomposition into plane waves (i.e. mimicking free photons) breaks down. Notwithstanding, the approximate formulae 
of Jones (1971) suffice for the purposes of this analysis (as will become evident shortly); the u 1 limit of his Equation (4) is 




do 

3T 


IB-WW 


7»1 , w«l , 


(18) 
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where we have included an extra factor of Z 2 to include consideration of ions other than protons. At photon energies 1 & u> 7 , 
the numerical computations in Haug (1972) indicate that the differential cross-section approximates a power-law (exceeding 
the analytic approximation of Jones 1971), before rolling over near the kinematic maximum energy of w ~ 7 . Equation (18) is 
readily integrated over the power-law Lorentz factor distribution for ions that can be obtained from Equation (15): 


dn^(u) 

dt 


inv 

brems 



r » ~ 1 
Tp — 1 + 3 tp 


pflV) 


,-r. 


U1 1 , 


(19) 


,(«) = f P log,(0.68i) = + j-L, . (20) 

This approximation to the contribution to the inverse bremsstrahlung emissivity from relativistic ions is still correct, up to a factor 
of order unity, when w ~ 1 , and hence the specification of the range w ^ 1 in Equation (19). When w ^ 1 , the spectrum breaks to 
assume a u>“( I+r p) form, so that the inverse bremsstrahlung spectrum is steeper than its bremsstrahlung counterpart by an index of 
roughly unity. We note that while the emissivity in Equation (19) can formally exceed that in Equation (11), in practice this never 
arises due to the distribution index T p exceeding 2 considerably at near-thermal energies, while T p £> 2 at relativistic energies. 
Hence, Equation (11) should be regarded as the appropriate form for optical and X-ray energies. A comparison of Equations (16) 
and (19) again yields the result that the factor [e p /e c ] (0T, P /0T,e) Tp ~ l roughly defines the ratio of inverse bremsstrahlung to 
bremsstrahlung emissivities at soft gamma-ray energies. 

It is also pertinent to briefly mention electron-electron bremsstrahlung. Since this process is a quadrupole interaction in a classical 
description, it is strongly suppressed for non-relativistic impact speeds relative to electron-ion bremsstrahlung. This is borne out 
in the differential cross-sections derived by Fedyushin (1952) and Garibyan (1953), which are of the order of da/du ; ^ 4ofro/(15o;) 
(see, for example, the exposition in the Appendix of Baring et al. 1999), a factor of j3 2 smaller than the result in Equation (4). 
However, the e — e bremsstrahlung cross-section for ultra-relativistic impact speeds is necessarily comparable to that for e — p 
bremsstrahlung (i.e. Equation [14]), as is indicated in the work of Baier, Fadin k Khoze (1967; see also the Appendix of Baring et 
al. 1999, and the discussion in Blumenthal k Gould 1970), who derived limiting forms using the Wei zs acker- Williams method. The 
similarity of cross-sections for 7 > 1 implies that the ratio of electron-electron bremsstrahlung to electron-proton bremsstrahlung 
is of the order of unity for plasmas with fairly normal ionic abundances, a fact that was noted by Baring et al, (1999). Hence, e — e 
bremsstrahlung can play an influential role in the determination of the importance of inverse bremsstrahlung. 

2.4. Expectations for Shock-Heated Environments 

Collecting all the rates so far assembled, it is now a simple enterprise to assert when inverse bremsstrahlung is significant in 
astrophysical scenarios. Setting 7£o-x to be the ratio of Equations (11) to (8) at w ~ (/? Tj e) 2 , so as to represent the ratio of 
inverse bremsstrahlung to bremsstrahlung at optical to X-ray energies, and #MeV to be the ratio of Equations (19) to (16) at 
w ~ 1 , representing the inverse bremsstrahlung/bremsstrahlung ratio at soft gamma-ray energies, we arrive at 


Hq-x 


£e 

Ce 


%) 


(Fp — 1)/2 


— ^MeV , 


(21) 


expressing the thermal speeds of the two species in terms of the electron and proton temperatures T e and T p . Similar order-of- 
magnitude estimates for such ratios can be obtained for other ionic species such as He 2+ . In Equation (21), it must be emphasized 
that for the optical/X-ray band ratio, the index T p refers to the mean index ranging from thermal proton speeds up to thermal 
electron ones , while for T^MeV , T p represents the mean index of protons with 7 1 . Hence, given the general upward curvature 

of spectra for cosmic rays accelerated in non-linear shocks (e.g. Eichler 1984; Jones k Ellison, 1991), we expect T p > 2 for the 
optical/X-ray band considerations and T p ^ 2 for the gamma-ray band. Since the gamma- ray spectrum for inverse bremsstrahlung 
(see the discussion after Equation [19]) is steeper than that for bremsstrahlung in Equation (16), 72-MeV defines an upper bound 
to the ratio of emissivities for the two processes in hard gamma-rays, a bound that is actually conservative given the expectation 
that e - e bremsstrahlung will amplify the bremsstrahlung signal by a factor of order 2 in the gamma-ray band. 

The ratio e p /e e of efficiency factors is a free parameter for a shock acceleration model. As noted at the beginning of Section 2, 
while efficient acceleration of protons is generally expected due to strong turbulence local to the shock (e.g. Ellison, Baring & 
Jones 1996), so that e p ~ 0.1 , to order of magnitude (which is borne out in observations of interplanetary shocks: see Baring et 
al. 1997), much less is known about the value of c e , essentially due to a paucity of in situ observations of electron acceleration in 
the heliosphere. However, this difficulty can be circumvented in part (at least in astrophysical scenarios) by arguing that cosmic 
ray data are a strong indicator of the acceleration efficiency in typical cosmic ray sources, principally supernova remnant shocks. 
The electron-to-proton abundance ratio above 1 GeV (where both species are relativistic) is well-known to be of the order of a 
few percent (e.g. see Muller et al. 1995), which constrains the typical acceleration efficiency, assuming that differences in electron 
and proton propagation in the interstellar magnetic field at these energies are not great. From the forms in Eq. (2), this ratio can 
be estimated as a function of the temperatures, masses and spectral indices of the two species. For an isothermal plasma with 
r c - T p ~ 2 , one quickly determines that n c (l GeV/c)/n P ( 1 GeV/c) - [e p /e e \ {m e /m p ) 1 ^ , so that the observed abundance ratio 
at momenta 1 GeV /c is reproduced only if e e is of the same order of magnitude as e p , i.e. injection of electrons into the Fermi 
mechanism is indeed efficient. Relaxing the isothermal and equal index stipulations yields a range of e p /c e on either side of unity. 
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As will be apparent in the subsequent subsection, a specific simulation of the acceleration mechanism applied to supernova remnants 
yields values of c p /e e not too disparate from unity. 

The temperature ratio T e /T p clearly controls the importance of inverse bremsstrahlung relative to bremsstrahlung, both through 
its explicit appearance in Eq. (21), and also via its influence on the efficiency ratio e p /c e . In shocked plasmas, this ratio is expected 
to far exceed m e /m p , thereby implying that inverse bremsstrahlung is insignificant in most astrophysical scenarios involving shocked 
plasmas . This is the principal result of this paper; it is contingent upon there being far fewer protons than electrons above electron 
thermal speeds, an occurrence we believe to be common near astrophysical shocks. 

Evidence for this can be derived from both observational and theoretical sources. On the observational side, there is a clear indi- 
cation of the lack of velocity equipartition in electron and proton populations near the Earth’s bow shock. Feldman (1985) exhibits 
measurements of electron and ion velocity distributions that lead to the identification of (i) strong electron heating during shock pas- 
sage, and (ii) downstream velocity thermal ratios for electrons to ions (mostly protons) of the order of 50 (i.e. m e T p /m p T e ~ 0.02 ). 
The situation for travelling interplanetary shocks is similar. While ion distributions indicate shock heating that couples to the 
ram pressure differential across such weak shocks (e.g. see Baring et al. 1997 for Ulysses SWICS data, where T p ~ 5 x 10 K), 
measurements of electron populations are generally sparse. Temperatures of the electron component of the solar wind upstream of 
interplanetary shocks is typically in the 3 x 10 4 -2 x 10 5 K range (Hoang et al. 1995), commensurate with large scale solar wind 
averages (Phillips, et al. 1995), and therefore are a sizeable fraction of the downstream ion velocity dispersion temperatures in 
the Ulysses data set. Hence, in the absence of any electron heating at the low Mach number interplanetary shocks encountered 
by Ulysses, m e T p /m p T e ~ 10~ 3 - 10~ 2 will be obtained. Such values are evinced in the ISEE data presented in Feldman (1985), 
which exhibits only modest heating of electrons by interplanetary shocks. 

Outside the heliosphere, the evidence for significant electron heating is sparse and more circumstantial. Shocks at the outer shells 
of young supernova remnants (SNRs) provide the best indication. These move typically at speeds 500-3000 km/sec, depending 
on the SNR age and the mass of the progenitor, suggesting ion heating to temperatures of the order of 5 x 10 7 -2 x 10‘ K per 
nucleon. Electron temperatures in SNRs can be deduced from observations of thermal X-ray emission (for a recent collection 
of observational studies, see Zimmermann, Triimper, & Yorke 1996), and are typically ~ 10 K. Therefore, it is expected that 
m e T p /m p T e - 5x 10“ 3 -0.2 in these sources. Inferences of electron to proton temperature ratios in active galactic nuclei are 
inconclusive, largely due to the absence of a marker for ion temperatures. 

The theoretical expectation that m e T p /m p T e < 1 is derived from considerations of dissipation in the shock layer. In a diffusive 
acceleration process where spatial diffusion is effected via elastic collisions of particles with magnetic irregularities that are anchored 
in the background fluid, the acceleration from thermal energies is effectively a velocity-dependent process (e.g. see the reviews of 
Drury 1983; Jones & Ellison 1991). This situation corresponds to dissipative heating in the shock layer that is independent of 
mass of the species. Baring et al. (1997) illustrates such a situation in their modelling of accelerated proton and He populations at 
nearby interplanetary shocks. However, dissipation in the shock layer cannot be entirely elastic in origin in any frame of reference, 
due in part to the different mobilities of electrons and ions. In a quasi-perpendicular shock, it is clear that the much smaller 
Larmor radii of electrons (relative to ions) that cross the shock lead to the establishment of induced electric fields via charge 
separation. In quasi-parallel shocks, similar electric fields exist, due to the turbulence of the field structure in the shock layer. In 
fact generalized shock structure including out-of-the plane magnetic field components requires the presence of electric fields in all 
frames of reference (Jones k Ellison 1987; 1991). Such fields heat the low energy electrons at the expense of the energy of the 
ions (see the discussion of Jones 1999), and thereby effect at least a partial equilibration in energy, providing significant theoretical 
support for the contention that m e T p /m p T e 1 . In general, a complex array of wave modes and instabilities can precipitate 
strong acceleration/heating of electrons. Limited evidence for such heating can be found in the hybrid plasma simulations of Cargill 
and Papadopoulos (1988) for non-relativistic shocks, and the full plasma simulations of Hoshino et al. (1992) for relativistic shocks 
in a pair-dominated ion-pair plasma. However, these two works do not simultaneously probe fully three-dimensional hydrogenic 
plasmas with truly realistic ion/electron mass ratios (so as to address a full array of wave modes). Three-dimensional simulations 
are required to accurately account for particle diffusion in all directions (e.g. see the discussion in Jones, Jokipii k Baring 1998). 
Hence, a definitive simulational assessment of the degree of electron heating in shocked plasmas remains an outstanding problem. 

2.5. A Simulational Illustration 

The expected dominance of bremsstrahlung can be illustrated using model predictions from a Monte Carlo simulation of diffusive 
acceleration in non-relativistic shocks. The Monte Carlo technique used to model diffusive shock acceleration is well-documented in 
the literature (e.g. see Jones k Ellison 1991; Baring, Ellison, k Jones 1993; Ellison, Baring, k Jones 1996; and more recently, Baring 
et al. 1999, for a description of electron injection). It is a kinematic model, closely following Bell’s (1978) approach to diffusive 
acceleration, where the simulation is used to calculate, in effect, solutions to a Boltzmann equation for particle transport involving 
a collision operator, without making any assumption concerning the isotropy of particle distributions. Particles are injected at a 
position far upstream and allowed to convect into the (infinite plane, steady-state) shock, diffusing between postulated scattering 
centers (presumably magnetic irregularities in the background plasma and self-generated turbulence) along the way. As particles 
diffuse between the upstream and downstream regions, they continually gain energy in accord with the Fermi mechanism. Particle 
splitting is used very effectively to maintain excellent statistics over a large dynamic range. The approach adopts a phenomenological 
mean free path, A, cc r| (with a - 1 for all ions and energetic electrons; for near-thermal electrons, see Baring et al. 1999 for 
details) to encompass the complications of plasma microphysics, where tv = pc/(ZeB) is the particle gyroradius. The aptness of 
this power-law prescription to shocked plasma environments is supported by particle observations at low energies at the Earth’s 
bow shock, Ellison, Mobius, k Paschmann 1990), deductions from ions accelerated in solar particle events (Mason, Gloeckler, 
&; Hovestadt 1983), and also from plasma simulations (Giacalone, Burgess k Schwartz 1992). Furthermore, in the strong field 
turbulence expected (and seen) at shocks, near-Bohm diffusion, A, £ r g , imposes an a — 1 requirement. 
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Fig. 1. — Accelerated particle distributions (a) from a Monte Carlo simulation of acceleration at steady-state cosmic-ray modified non-relativistic 
shocks, specifically from the modelling by Ellison, Goret & Baring (1999) of broad-band spectra for the supernova remnant Cas A, together with (b) 
the bremsstrahlung and inverse bremsstrahlung emissivities resulting from these distributions. The electron and proton distributions are omni-directional 
fluxes, binned in kinetic energy intervals, in accord with our previous expositions, and represent the differential energy distributions multiplied by the 
particle speed; they were obtained for a shock speed of 2800 km/sec and an ISM density of 3cm‘ 3 . An approximate power-law velocity distribution with 
index T p = 3.6 is depicted to facilitate the discussion in the text. The photon emissivities are multiplied by photon energy, i.e., E y times the differential 
spectrum. 

Typical predictions for the proton and electron distributions (downstream of a shock) using the Monte Carlo technique are 
depicted in Figure 1, taken from the Ellison, Goret & Baring’s (1999) study of shock acceleration in the supernova remnant Cas 
A (similar spectra are presented in Baring et al. 1999). The distributions are exhibited as omni-directional fluxes (i.e, differential 
energy fluxes multiplied by particle velocity: vdn/dE = c0n(E )), mimicking the spectra measured by particle detectors in the 
heliosphere. Hence, for non-relativistic speeds, the momentum (or velocity) distribution index is exactly twice the slope of 
the curves in Figure la, while for ultra-relativistic particles, T s is identical to the depicted slopes, which for both the electrons 
and protons is roughly T e « T p « 2.2 . The particular simulation output that is illustrated possesses a hot electron component 
in the downstream plasma, with a temperature roughly half that of the protons. The light power-law drawn in Figure la, with a 
velocity distribution as labelled, approximates the proton distribution between thermal proton and thermal electron speeds; it has 
an index of T p « 3.6 (i.e. flux distribution 0 p n p {E p ) oc E~ l s ). Observe that the Monte Carlo curves are not pure power-laws in 
momentum but exhibit a slight upward curvature in the relativistic portions of the distributions (and also significant curvature at 
suprathermal non-relativistic energies for the protons), a central feature of the non-linear feedback between the ion diffusion and 
the shock hydrodynamics in cosmic-ray modified shocks that is well-documented (e.g. Eichler 1984; Ellison & Eichler 1984; Jones 
& Ellison 1991; Ellison, Jones k Baring 1996). Because such hydrodynamics is strongly regulated by the feedback between particle 
acceleration and escape from the remnant, the generated proton distributions are not very sensitive to the prescription of A as a 
function of rigidity. 

Electrons, on the other hand, are dynamically unimportant, so that their spectral normalization (but not spectral shape) is 
sensitive only to their injection efficiency, which couples to the electron-to-proton temperature ratio. While this is somewhat 
uncertain, we used reasonable shock speeds and observed X-ray temperatures from remnants (e.g. Zimmermann, Triimper, k 
Yorke 1996) to adopt appropriate values. Note also that the ratio of the electron to proton populations in the 1-10 GeV range 
(contingent upon the assumed injection efficiency) for the specific results in Figure la varies between 15% and 5%, somewhat larger 
than observed e/p cosmic ray ratio in this energy range (e.g. see Muller et al. 1995); e/p ratios more consistent with the cosmic 
ray data are readily obtained (Baring et al. 1999) by the Monte Carlo simulation. 

Figure lb depicts the emissivities for bremsstrahlung (solid curve) and inverse bremsstrahlung (dashed curve) that result from 
the particle distributions in Figure la. The flat spectra below X-ray energies correspond to dn 7 (u/)/<ft cc w” 1 , approximating the 
low energy results in Equations (8) and (11), and interestingly, close to the index observed in the cosmic X-ray background below 
10 keV (thus motivating the application discussed in Boldt k Serlemitsos 1969). Above the electron thermal energy, the spectra 
steepen for both processes, and in the gamma-ray regime, the inverse bremsstrahlung spectrum (oc 2 « ) is steeper 

than that of bremsstrahlung ( a a? -2,2 « w ~ T<l ), in accord with Equation (16) and the discussion following Equation (19). 

We note that the steep shoulder in the bremsstrahlung spectrum in the X-ray band defining the roll-off above the thermal emission 
could potentially provide an alternative explanation (Baring et al. 1999) to synchrotron emission for the non-thermal X-rays seen in 
the supernova remnants SN 1006 (Koyama et al. 1995), IC 443 (Keohane et al. 1997) and Cas A (Allen et al. 1997). The illustrated 
case, of index ~ 2.8 for just over half a decade in energy, lies in between the index of ~ 2.3 ASCA obtained for IC 443 and the 
X-ray observations of Cas A and SN1006 (~ 3). While this issue is discussed in more detail in Ellison, Goret k Baring (1999), it 
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is important to remark that inferences of the spectral index of a non-thermal component in X-ray data are very sensitive to the 
choice of temperature of a thermal component and the shape of a non-thermal continuum. Restrictions to pure power-laws for 
the continuum, while convenient, are not really appropriate for either synchrotron radiation turnovers or bremsstrahlung emission 
mechanisms, so that conclusions about the physical mechanism in operation should be guarded. For example, the RXTE data 
in Allen et al. (1997) on Cas A suggests a spectral break at 16keV, steepening from an index ~ 1.8 to one of - 3.0, data 
that at first sight does not appear inconsistent with a bremsstrahlung model in the transition region from thermal to non-thermal 
energies (though probably with a slightly higher temperature than the case depicted in Figure 1). The critical distinction between 
synchrotron turnover and “transition” bremsstrahlung models therefore lies in the expected flattening above the transition region 
in the bremsstrahlung spectrum (e.g. see Figure lb). This provides an observational discriminant via broader-band observations 
beyond the scope of ASCA and RXTE, though current instrumental sensitivities (such as those of XMM and Integral) are probably 
insufficient to answer this question. A related situation arises with diffuse/unresolved emission in the galactic ridge, and recent 
broad-band X-ray spectra of the Scutum arm (Valinia, Kinzer & Marshall 1999) suggest a turnover in the 50-100 keV range, which 
could indicate the superposition of Comptonized (or synchrotron) emission from unresolved discrete sources. However, we note 
that data statistics in the OSSE band are not sufficient to preclude a flattening at higher energies that would be characteristic of 
bremsstrahlung of either discrete source or diffuse origin. 

Returning to the issue at hand, the particular example in Figure 1 gives a strong dominance of bremsstrahlung over inverse 
bremsstrahlung. In fact, given the r p « 3.6 index for non-relativistic speeds, and T e /T p ~ 1/2 , Equation (21) can be simply 
applied to estimate a ratio of ~ 8000 between the two processes (for e p /€ e = 1), close to the ratio deduced in the optical band 
from Figure lb. This agreement highlights the viability and usefulness of the schematic distributions in Eq. (2), and also that of the 
consequent quantitative estimates of the various emissivities throughout this paper. Observe that the only way to increase 77o-x 
is to dramatically lower the temperature of the electrons relative to the protons, but to do so would cause a corresponding decrease 
in the e/p ratio at relativistic energies (and influence the ratio e p /e e of injection efficiencies), well below observed values in the 10 
GeV band (see, e.g. Muller et al. 1995, for a recent compendium of data). As argued above, it is unlikely that fto-x ^ 0.01 will 
be realized often in shocked astrophysical plasmas. 

It is natural to question whether such a result holds for ions heavier than protons. Clearly, both the bremsstrahlung and inverse 
bremsstrahlung contributions from helium and heavier elements are given by their abundance (usually small) multiplied by Z 2 . 
The Z 2 enhancement provides significant bremsstrahlung emissivities when helium targets (typically of ~ 5% abundance) are 
involved. For inverse bremsstrahlung, the Z 2 amplification potentially could be offset by a more severe mass ratio reduction, if 
lower thermal speeds are precipitated downstream of shocks. There is little or no evidence for such lower thermal speeds for ions 
heavier than protons. In fact, the Ulysses SWICS data on the 91097 interplanetary shock unequivocally indicates the contrary: the 
shock dissipation mechanism generates virtually identical velocity dispersions for both protons and alpha particles (Baring, et al. 
1997). Such a property is well-modelled by the Monte Carlo technique, since the elastic scattering hypothesis that it employs renders 
the shock acceleration model a velocity-dependent injector. This scattering assumption, together with a rigidity-dependent diffusion 
coefficient, is rigorously put to the test in the successful prediction (Meyer, Drury, & Ellison 1997, Ellison, Drury, & Meyer 199/) 
of abundances of ions of various metallicities in the cosmic ray population, using a supernova remnant site for the acceleration of 
dust grains. These pieces of evidence point to an approximate velocity equipartition of thermal ions in shocks, so that the criterion 
in Equation (21) renders the conclusion that inverse bremsstrahlung is generally insignificant in shocked astrophysical plasmas also 
applicable to contributions from ions heavier than protons. 

3. DISCUSSION 

The unimportance of inverse bremsstrahlung for shocked environs in discrete sources, the principal result of this paper, can 
be interpreted in the broader context of invocations of inverse bremsstrahlung in diffuse emission problems. The pre-eminent 
implication of this result is obviously that supernova remnants will dominate their environs in X-ray and gamma-ray luminosity, 
and therefore will obscure truly diffuse emission components in extended regions such as the galactic ridge and the Orion complex, 
if the source remnants are too densely clustered. Hence, it is of significant interest to identify the length scales on which discrete 
sources such as supernova remnants provide a dominant contribution to unresolved X-ray and gamma-ray emission in the galaxy, 
as a guide to modelers and also to aid the interpretation of experimental results. 

The discussion of this issue is not restricted to just the ratio of bremsstrahlung from cosmic ray primary (i.e. shock-accelerated) 
electrons and inverse bremsstrahlung from cosmic ray ions. The works of Tatischeff, Ramaty &: Kozlovsky (1998) and Dogiel et al. 
(1998) argued that suprathermal proton bremsstrahlung emission levels are comparable to those of bremsstrahlung from knock-on 
electrons in the Orion complex, given the then (and now defunct) Comptel nuclear line detection in this region. Knock-on electrons 
are generated from the cool ambient interstellar population via Coulomb collisions with energetic cosmic ray ions (and perhaps 
electrons). Since the stopping length (on ambient electrons or protons) for 25 keV electrons in interstellar environments of typical 
density is of the order of 100 pc (e.g. see Valinia & Marshall 1998), Coulomb collisions are clearly relevant for electrons in large 
scale systems such as the Orion cloud complex. Coulomb scattering simultaneously mediates both energization of ambient electrons 
and cooling of cosmic ray primary electrons. In contrast, shocked astrophysical plasmas such as those in young to middle-aged 
supernova remnants are collisionless, so that knock-on populations are locally negligible. Given the conclusions of this paper 
relating to shock-heated plasmas, it is salient to also examine when knock-on electron bremsstrahlung is a significant contributor 
to diffuse/unresolved emission relative to primary electron bremsstrahlung. Note that the comparisons of Tatischeff, Ramaty & 
Kozlovsky (1998) and Dogiel et al. (1998) of inverse bremsstrahlung relative to knock-on electron bremsstrahlung did not include 
an abundant, shock- accelerated electron population, and therefore did not address this issue. 
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3.1. Knock-on Electrons 

First, we address the question of whether knock-on electrons are abundant relative to cosmic ray primary electrons. For mono- 
energetic cosmic rays, non-relativistic knock-on electrons possess an E~ 2 kinetic energy distribution regardless of whether the 
energetic cosmic rays are ions (e.g. see Hayakawa 1969) or electrons (Gould 1972; true also for relativistic electrons: see Baring 1991), 
which corresponds to a /?” 3 velocity distribution. This is comparable in steepness to the shock- accelerated electron distribution 
exhibited in Fig. la. However, in spite of the ambient (and therefore also knock-on) electrons being much more abundant than 
cosmic ray ones, the relative population level at almost relativistic speeds defines the criterion for the importance of knock-on 
electrons in producing X-ray bremsstrahlung. The knock-on electron population extends upward from extremely low thermal 
energies (temperatures of 3-100 K; see Dogiel et al. 1998 for a listing of parameters for the Orion complex). Assuming that only 
a minority of ambient electrons can be energized in knock-on collisions without disturbing the thermal balance of the ambient 
interstellar medium, it is straightforward to determine that the distribution n e (£ , Mev) ~ 10' -»n.TiEj£ v MeV 1 cm 3 forms a 
conservative upper bound to the knock-on energy distribution out to MeV energies. Here 7\ is the ISM temperature in units of 
10 K, and n e is the ISM density in units of cm’ 3 , typically of the order of 1 cm" 3 in the normal ISM, and ~ 100 cm" 3 in dense 
interstellar regions. Also, Em*v is the electron energy in units of MeV. 



Fig. 2. — A schematic depiction of the knock-on electron distribution (heavy line) generated by a cosmic ray ion distribution that is a power-law above 
thermal energies represented by a non-relativistic temperature 7\ ( <£ m.c 2 /^). For ions originating in SNRs, typically kTi ** 10“ - 1 MeV from shock 
heating. The knock-on distribution extends upwards from the temperature Tism of the ambient ISM, and steepens above the kinematic break at twice 
the cosmic ray ion thermal speed (i.e. E ** (4‘me/T7i^)fcT , » j see Baring 1991, for example, for a discussion of kinematics). While knock-on electrons may 
dominate cosmic ray primaries by number, they suffer a substantial paucity at high energies. 

This upper bound is actually a gross overestimate, given that it corresponds to 100% of the ambient electron population being 
energized. Furthermore, due to the kinematic maxima for the energies imparted to cold electrons by Coulomb collisions, above the 
thermal speed of the ballistic cosmic rays, which for ionic projectiles is unlikely to exceed typical shock speeds for supernova remnants 
of a 500-3000 km/sec, the knock-on spectrum steepens from the dependence E ~ 2 . The knock-on distribution can then be shown 
to roughly follow an £“ 3 / 2_r p distribution for a cosmic ray ion spectrum a E~ r * (derivable, for example, using Eq. (9) of Baring 
1991), rendering the knock-on electron population at around 100 keV far smaller than the above bound; see Figure 2 for a schematic 
depiction. An approximate refined estimate for the knock-on distribution would then be n e (FJ M ev) ~ 0.1n e T\ [E’Mev/lO -4 ] -3 ^ 2- ^ 
MeV -1 cm -3 for £?MeV 10“ 4 . This can be compared with a rough estimate for the non-thermal electron distribution expected 
for efficient acceleration at a remnant’s shock: n e (£ , M ev) ~ 10 5 n e T 7 [E M ev/10' 4 ]" 2 MeV- 1 cm” 3 , where here T 7 is the shocked 
plasma temperature in units of 10 7 K. Since X 7 ~ 1 near SNR shocks, clearly the primary cosmic ray source spectrum locally 
swamps the knock-on one. Away from such shocks, diffusion dilutes the primary source spectrum down to cosmic ray flux levels. 

While a specific prediction for the ambient cosmic ray primary electron density could be made via a particular acceleration model 
(such as the IMonte Carlo technique used above), combined with a specific diffusive propagation model, we prefer to make use of 
existing experimental cosmic ray data under the premise that dense interstellar regions are not extraordinarily peculiar in their cosmic 
ray properties. Muller et al. (1995; see also Muller & Tang 1987) exhibit the primary electron flux at Earth above around 3 GeV. 
From their data, one quickly arrives at the differential energy distribution for electrons of n fi (£ M ev) ~ 10 -4 #MeV 5 MeV-1 cm_3 afc 
energies above 3 GeV. Direct extrapolation of this down to 100 keV (i.e. well-below energies where heliospheric modulation impacts 
measurements of the cosmic ray spectrum) would lead to the conclusion that cosmic ray electrons should dominate the maximum 
possible density of knock-on electrons by around 4 orders of magnitude, and the refined estimate above of the knock-on density by 
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a factor of 10 9 or more (for T p £ 2), for realistic values of the ISM temperature T x . However, the primary electron spectrum is 
probably much flatter in the 100 keV to 3 GeV band than measured at higher energies. The diffuse radio synchrotron emission in 
the 30-600 MHz band (a v ~ 2/3 ) in our galaxy suggests (e.g. see the review by Webber 1997) a cosmic ray electron distribution 
<x E~ 7/3 between 500 MeV and 2.2 GeV. This flattening probably reflects the transition from directly sampling contributions to 
the cosmic ray spectrum near discrete sources (which generally should possess flatter — E~ 2 distributions like those in Figure 1) to 
larger-scale steepening imposed by cosmic ray propagation effects. Even using an extrapolation n e (i?Mev) ~ 10“ E MeV MeV" 
cm" 3 at energies below 300 MeV based on the inferences from the radio data, one arrives at the conclusion that the cosmic ray 
primary electron population would far exceed any knock-on component generated by cosmic ray ions. 

3.2. Diffuse and Unresolved Emission 

Hence, it follows that a knock-on population of electrons can only generate a significant amount of electron bremsstrahlung 
emission in dense interstellar regions provided that there is a profound paucity of cosmic ray electrons in the 30 keV-1 MeV range. 
Cooling of electrons through Coulomb and ionization losses potentially can create such a paucity; lower energy electrons cool 
faster than more energetic ones do. Therefore, a depletion of cosmic ray electrons below 1 MeV is a distinct possibility on scales 
considerably greater than lOpc outside the supernova remnants that act as sources of cosmic rays. Specifically, using the numbers 
quoted by Valinia and Marshall (1998), 25 keV electrons that would emit bremsstrahlung in the RXTE band have a collisional loss 
pathlength of around lOOpc in the ISM. This linear scale A Cou ! can be translated into a radial scale r Co u! by incorporating the 
effects of diffusion due to the turbulent interstellar field. The mean free path A B for such diffusion is extremely uncertain, but as is 
the case for the undisturbed interplanetary medium in the heliosphere, estimates for A B (by coupling to the field irregularity scale; 
e.g. see Kim, Kronberg k Landecker 1988) far exceed particle gyroradii, at around lOOpc. This implies that for 25 keV electrons, 
rcoui = [Acoul A b ] 1/2 ~ lOOpc defines the maximal extent of their “Coulomb halo” around a SNR. For higher energy electrons, 
the radius of this halo is larger, so that at any given radius, the primary electron spectrum is depleted up to the energy at which 
Coulomb losses become insignificant. 

Irrespective of whether and on what scales such depletions exist (no data pertaining to cosmic ray electrons samples this energy 
range), the dominance of bremsstrahlung from primary cosmic ray electrons over inverse bremsstrahlung from primary cosmic ray 
protons that is established in this paper would lead to the major contribution to unresolved “diffuse” X-ray emission in Orion or the 
galactic ridge coming from discrete sources, unless the volume filling factor of supernova remnants in the interstellar medium is less 
than 10" 4 -10“ 3 . While accurate determination of this factor is difficult, extant compendia of supernova remnants (e.g. Helfand 
et al. 1989; Green, 1999) lead to estimates at the lower end of this range (i.e. ~ 10“ 4 ), to order of magnitude. Hence, while it 
is possible that such conditions might prevail in some galactic interstellar environments, and more probably at moderate latitudes 
above the disk, we contend that primary cosmic ray electron bremsstrahlung will simultaneously dominate bremsstrahlung from 
secondary knock-on electrons and inverse bremsstrahlung mediated by cosmic ray ions as contributors to diffuse X-ray emission in 
the most important and interesting environments. 

When this is not the case, primary electron bremsstrahlung signals (of spectral index 2-3) will be confined to discrete source 
locales of scales ~ 10 -lOOpc, i.e. near supernova remnants, while suprathermal proton bremsstrahlung and knock-on electron 
bremsstrahlung (of similar spectral indices) will sample the length scales appropriate to ionic collisional losses, namely ~ 0.2 Mpc 
for the 50MeV protons that can generate inverse bremsstrahlung X-rays in the RXTE band. To reiterate, these length scales will 
be appropriately contracted in linear dimensions if they exceed the diffusive scales for such particles in the interstellar magnetic 
field (e.g. to - 5kpc from the protons: see Valinia and Marshall 1998). Such radiation “clumping” issues should form a focus 
of discussions of diffuse and unresolved emissions. Note that because - 50 MeV cosmic ray protons are not easily depleted on 
scales 5kpc, and since the low ISM temperature forces the knock-on distribution to low levels, we expect that generally inverse 
bremsstrahlung dominates knock-on electron bremsstrahlung outside remnants. 

An observationally interesting question concerns how easily bremsstrahlung signals from primary cosmic ray electrons associated 
with SNRs can by distinguished from truly compact sources such as neutron stars or black holes, given current angular resolution 
capabilities of X-ray telescopes. Spectral signatures may provide the answer prior to improvements in instrumental resolving 
power. The similarity of unresolved emission spectra in source-rich and source-poor galactic regions (ASCA data: Tanaka, Miyaji 
& Hasinger 1999), and the spectral similitude of observations of the Scutum arm (RXTE/OSSE data: Valinia, Kinzer k Marshall 
1999) and galactic X-ray binaries and black hole candidates provide clues impinging upon this issue. 

Observe also that the expected general dominance of bremsstrahlung from primary electrons returns one to the energetics issues 
raised by Skibo, Ramaty k Purcell (1996), who discerned that the general inefficiency of bremsstrahlung processes in the X-ray 
band relative to collisional and ionization losses elicits difficulties in generating the requisite luminosity. In passing, we note that 
clumping of emission in the proximity of remnants will help alleviate the energy supply problem for a bremsstrahlung origin of the 
galactic diffuse X-ray emission, by diminishing such losses for the electrons. 

4. CONCLUSION 

In conclusion, the results of this paper indicate that inverse bremsstrahlung in the optical to X-ray and gamma-ray bands can be 
safely neglected in most models of discrete sources invoking shock acceleration. This follows from the approximate ratio, posited 
in Equation (21), of inverse bremsstrahlung to bremsstrahlung from primary non-thermal electrons, together with expectations for 
dissipational heating of electrons in the shock layer. Both observational and theoretical evidence favors electron temperatures almost 
comparable to, and certainly not very deficient relative to, proton temperatures in shocked plasmas. Hence shocked environments 
in discrete sources such as supernova remnants should have inverse bremsstrahlung contributing only in a minor capacity to the 
overall emission, a situation that could apply also to larger, extended interstellar regions such as the galactic ridge, provided that 
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the volume filling factor associated with remnants exceeds around 1O~ 4 -10" 3 . Since this filling factor may sometimes be lower, 
perhaps at moderate to high galactic latitudes, we find that in such cases, inverse bremsstrahlung can contribute significantly in the 
X-ray band on length scales of the order of 5O-10Opc or greater, thereby posing an interesting observational goal to resolve bright 
bremsstrahlung X-ray emission concentrated around discrete sources from any truly diffuse inverse bremsstrahlung component that 
spatially traces the cosmic ray ion population. We also contend that bremsstrahlung from knock-on electrons is probably always a 
minor contribution, and find that the dominance of bremsstrahlung from primary cosmic ray electrons over inverse bremsstrahlung 
in the vicinity of remnants is even more enhanced in the gamma-ray band. 

We thank Stephen Reynolds, Elihu Boldt, Azita Valiniaand Peter Serlemitsos for discussions, Vincent Tatischeff and the anony- 
mous referee for comments helpful to the improvement of the presentation, and Bob Gould for pointing out certain literature on 
bremsstrahlung to us. 


APPENDIX 


In this brief Appendix, a pedagogical alternative to the derivation of the various bremsstrahlung emissivities is provided, using 
approximate, simplified cross-sections, so as to elucidate the results obtained in the text by diminishing the mathematical complexity. 
To this end, the expressions obtained here are proportionalities, possessing only the principal variables. Eqs. (1) and (2) are again 
the starting points. Instead of using Eqs. (4) and (14) for the non-relativistic and ultrarelativistic bremsstrahling differential 
cross-sections, it is expedient to use a single approximation 


da 

did 


NR-ER 



2 r n 



(Al) 


where 0(z) is a step function that is unity for 0 < x < 1 , and zero otherwise. Consider first normal bremsstrahlung. The 
integrations over the electron distribution are straightforward, with both the p < p T ,e and p > p r>e portions contributing for 
u> e <£ 1 , and just the high momentum tail being sampled for «w < 1. The resulting proportionalities can then be 

summarized as (for A = Z 2 n p n e oc{r^c) 


dn 7 (u) 

dt 


brems 


A 

0T,e 
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(A2) 


dependences that are clearly evident in Eq. (8); the absence of logarithmic factors is an obvious consequence of the simplification 
of the cross-section. There is little need to explicitly state the equivalent proportionality for the case of non-relativistic inverse 
bremsstrahlung, which reproduces the principal components of Eq. (11) in a similar manner. For normal bremsstrahlung from 
ultrarelativistic electrons, the integration can be performed with equal efficacy using Eq. (15), yielding 


dn 7 (u;) 


dt 


(0r, 


Ae« 
T=r e 


-r e 


1 £ 


(A3) 


thereby reproducing the form of Eq. (16), minus the complexity of integrations over the combinations of logarithmic and power-law 
factors. 
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